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Abstract 

Ballistic transport of electrons through a quantum wire with a constriction is 
studied in terms of Bohm's interpretation of quantum mechanics, in which the 
concept of a particle orbit is permitted. The classical bouncing ball trajecto- 
ries, which justify the name "ballistic transport" , are established in the large 
wave number limit. The formation and the vital role of quantum vortices is 
investigated. 
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The ballistic transport of electrons in quantum devices based on the two-dimensional 
electron gas is a subject of intense study. The name "ballistic transport" is borrowed from 
classical mechanics where a point-like particle moves freely in a region except for elastic 
collisions with obstacles or the confining walls. Obviously, the classical orbits are straight line 
segments connected through a reflection law. Quantum devices based on ballistic transport 
can be approximately thought of as an electron waveguide, in which the electrons behave 
as non-interacting particles of effective mass m* in a constant potential with hard wall 
confinement. While the picture of bouncing balls is a vivid one, quantum mechanically we 
need a precise justification for this picture. For this purpose, we study here a very simple 
model, that of a quantum wire with a constriction, to see how classical orbits emerge from 
the quantum wavefunction. 

We consider a straight quantum wire of width d. At the origin x = we place a thin 
but opaque obstacle attached to one side wall so that the channel opening is narrowed to a. 
Electrons of energy E are injected from the left in transverse mode n,. In order to develop 
our argument, we solve this problem by the mode matching method. 

Denote the wavefunction to the left and right of the narrow constriction as and 
respectively. These can be expanded in plane waves as 

*r = E^e^ + C Tn e- ianX ) sm(k n y) (1) 

n 

where r = 1, 2, k n — nn/d, and 



a., 



= ^E-kl (2) 



with appropriate units such that h 2 /(2m*) = 1. We also define symbols Cp n = Cm ± Cm- 
At x = 0, the wavefunction in the opening is expanded as 



sn(y) = V\ x=0 



£ mJ D m sin^ 0<y<a 

(3) 

a <y < d 

Continuity of the wavefunction at x = requires that 



C\ n — _ X! InmD m (4) 

rra 

where 



2 r a wry miry 

Inm = -i sin — — sin dy 

a Jo a a 

2 Tsinr^f - m)a sin(*f + m)a 



mr m7r nn i m,7r 

da d a 



(5) 



Similarly, continuity of the wavefunction derivative at x = results in 

(6) 

ra ra 

In matrix form, Eqs. (|4]) and (|6]) read: 

C+ = C+ = ID (7) 
/*aC{- = PaC 2 (8) 

Notice that the above equation does not imply Cf = C%, since J is not invertible. 

Physically, the scattering boundary conditions require that C 2 = so that no backward 
wave is traveling in region 2; therefore 

C 2 = ID and C x = ID - C x . (9) 

Inserting this into Eq. @, we obtain 

D = (Pal^Pad (10) 



-'raraj 



Thus the wavefunction is determined by the injection mode expansion coefficients C\ n = 5 n 

Calculating the transmission coefficient is not the main purpose of this paper. Rather, 
we wish to examine how the electron passes through the channel, and to compare this 
with the corresponding classical motion. To this end, we adopt the formalism of Bohm's 
interpretation of quantum mechanics HUH, in which a system is described not only by its 
wavefunction, but also particle orbits. In this picture, a particle at position r moves with 
velocity v = j/|\l/| 2 , where j is the probability current density. Since Bohm's picture is most 



vividly applied to scattering states, Holland in his recent book called it "the quantum theory 
of motion" [[J. We are dealing with a time independent wave equation, so the process of 
finding the particle trajectories is separated from solving for the wavefunction. All we have 
to do is to find the stream lines of the velocity field deduced from the wavefunction; they 
are curves tangent to j. 

According to our recent discussion of quantum probability flow 0|| , the stream lines can 
be classified in three categories: 1) those that close on themselves; 2) those which participate 
in global transmission, (which start and end at infinity; they could be regarded as special 
stream lines of type 1 which close at infinity;) 3) those which connect saddle points or end 
on the boundary; these lines separate different pieces of the flow. In order to study particle 
transport we concentrate on type 2 stream lines as shown in Fig. 1. We drew these lines 
by starting from points along the channel opening at x = 0, and tracing the trajectories 
forward and backward. The density of the starting points was taken to be proportional to 
the probability density in the opening. In order to exhibit the relation to the classical orbits, 
we deliberately chose a high electron energy. The trajectory in the transmitted region, as 
expected, clearly resembles that of a classical ball bouncing along between hard walls, or light 
rays between parallel mirrors. We can understand this better by developing an approximate 
analytical solution. 

Consider a situation where the transverse and longitudinal wave numbers are all large. 
The computed transmission coefficient is very close to the classical prediction, a/d. A 
large transverse wave number implies that I nm in Eq. will have appreciable non-zero 
values only when mr/d pa mir/a. From this condition and Eq. (|T0|), we see that D m is 
concentrated in a narrow interval centering about m = rii(a/d). Also from Eq. (|9|), C^n and 
Ci n are concentrated about n = rii. These estimates are well verified by exact numerical 
calculations. Denoting k x = a ni , k y = k no a n can be expanded about k x : 

«n - k x = -(k n - k y ) tan# (11) 

where 
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tan# = k y /k x . 



(12) 



Here we have used the condition of large longitudinal wave number. Under this approxima- 
tion, 

\I>2 ~ exp[i(k x + k y tan6)x] x 

Cm exp(—ik n x tan 6) sm{k n y) 

n 

= exp[i(a ni + k ni tan#)x] x 

^ cos k n (y — xta.n9) — cos k n (y + xtani 



2i. 



sin k n (y — xtan6*) + sin k n (y + xtan6>) 



(13) 



Expressing cos k n (y ± x tan#) as (l/k n )d[sm k n (y ± x ta.ia9)}/dy, and approximating l/k ri 
1/ky in the summation, we obtain 

^2 — exp[i{k x + kytan6)x] x 
f 1 d 



2ik y dy 



sn(y — x tan 9) — sn(y + x tan ( 



+ - [sn(y — x tan 6) + sn(y + x tan 6)] j 



(14) 



where the function sn(y) is defined by Eq. @, together with its odd symmetry extension, 
and its periodic extension with period 2d. In other words, 



sn(y) 



E m D m sin 



in-y 



\y\ < a 



(15) 



\y\ > a 

where y is y minus an integer multiple of 2d such that the result falls in (— d, d]. Now 



return to Eq. (|14D and consider the derivative of the function sn(y). Its non-zero part is 
J2m D m (mir/a) cos[m7ry/a\. Now since m is in a narrow region of riia/d, the factor mn/a 
may be approximated by k y (this kind of approximation would not be permitted when the 
factor is inside a rapid changing function such as sin, cos, and exponential). This leads to 
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\& 2 ~ exp[i(£; x . + k y tan8)x] x 

— [cs(y — x tan 9) — cs(y + x tan 0) 
+ i sn(y — x tan 9) + i sn(y + x tan 9)} 

= exp[i(k x + k y tan9)x] x 

ep[z"(y — xtan#)] — ep[— i(y + xtan( 
2i 



(16) 



where 



cs(y) 



EmDmCOS inm \y\ <a 



\y\ > a 



(17) 



and 

ep(iy) = cs(y) + isn(y) . (18) 

Aside from a phase factor, Eq.(^) expresses ^2 i n terms of an exponential-like function 
ep(iy). When x = 0, it clearly reproduces Eq. (|^) since sn(y) = [ep(iy) — ep(—iy)]/(2i). 
As x increases, ep[i(y — x tan 9)} represents a function ep[iy] shifted by xtan# upwards, and 
ep[—i(y + xtan#)] represents a function ep[— iy] shifted the same amount downwards. In 
other words, the two functions ep(±iy) each move at an angle ±9 to the channel walls as x 
increases. The motion of the functions leaves stripes of non-zero parts of the wavefunction. 
Fig. 2 illustrates the formation of the bouncing ball pattern in between the two thick 
horizontal lines which represent the wire walls. At the left edge of the figure, vertical line 
segments locate the non-zero parts of the function ep(±iy). The motion of these vertical lines 
traces out the upward and downward inclined stripes. What has physical meaning are those 
portions where the stripes lie inside the wire. Regions marked by a right upward arrow (U- 
regions) are those due to ep(iy), and those marked with a right downward arrow (D- regions) 
are due to ep(— iy). The triangular regions marked with horizontal arrows (H-regions) are 
where the two contributions superpose. 
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Not only can we show that the non-zero part of the wavefunction produces the pattern 
shown in Fig. 2, we can also show that the trajectories indeed follow the directions shown. 
To see this we calculate the particle velocity 

v = (h/m*)lm(^*V .^)/|^| 2 (19) 

To compute the gradient of ^2, first note that 

^M=i*„ep(%) (20) 

in the approximation mn/a ps k y that we have used before. Thus ep(iy) behaves much 
like exp(ik y y) in the non-vanishing region. Consider a U-region, where \&2 = exp[i(k x + 
k y tan 9)x]ep[i(y — xtan9)]/(2i), 

' 1***2, = ik y * 2 (21) 



X - 4 ? 

ox ay 



which lead to 



v x = (h/m*)k x , v y = (h/m*)k y ; (22) 



therefore 



v y /v x = tan 9 . (23) 

That is to say the particle moves to the right and upward in a direction determined by the 
decomposition of the total wave number in the x and y directions. Similarly the particle in 
a D-region moves in the direction v y /v x = — tan 9. 

In an H-region, it is easily seen that v x takes the same value as in the previous two cases. 
For v y , one finds 

ik y /4:{ep*[i(y — :rtan#)] — ep*[— i(y + xtan9)]} x 

{ep[i(y — xt&n9)} + ep[— i(y + xtan9)]} (24) 



which yields 

HO- 

k y /4 [\ep[i(y - x tan 9)]\ 2 - |ep[% + :rtan6>)]| 2 } . (25) 

Eq. (2"5|) is a small quantity due to the property of the exponential-like function ep(iy). 
The H-regions are symmetrical triangles. Obviously, v y vanishes when close to the bottom 
side of the triangle. Along the symmetry line of the triangle, xtaxiO = (integer) x d, which 
implies |ep[i(y ±xtan#)]| are equal; thus v y = 0. At other points of these triangles, v y takes 
on small values, which allow a smooth connection of the trajectories to the other regions. 

So far we have shown clearly the formation of the bouncing ball trajectories in the 
downstream transmitted region. Let's now discuss the reflected region. From Eq. ([J), (|9D 
and Cin = 5 nni , one has 

*! = e ik * x sm(k n y) + £(C 2n - C ln )e~^ x sm(k n y) 

n 

= 2i sm(k x x) sm(k n y) + ^ 2 {~x, y) 

= ^ R + ^ 2 (-x,y) (26) 

Here is a standing wave solution for the total reflection problem as if the channel were 
completely closed at x — 0; it represents a standing wave in the upstream or left part of the 
channel. Thus the complete wave function upstream is a superposition of this standing wave 
and a reflection of the transmitted wave. The standing wave has a grid of nodal lines which 
divides the region into cells of size Ax = n/k x , Ay = ir/k y . The magnitude of the standing 
wave at the center of each cell is ±2. This pattern will be preserved in at the places not 
covered by ^2{—x,y)- The velocity in these places, however, will be virtually zero since 
is pure imaginary. 

At the places covered by ^(^ x,y), the interference with the standing wave can form 
quantum vortices [[|. These quantum vortices perturb the stream lines, causing many 
detours away from straight line trajectories. This is why the trajectories upstream are 
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much more complex than those downstream in Fig. 1. Quantum vortices appear in many 
electron waveguide structures, for examples, in that of Lent , and Berggren |7| . 

We now explain the deployment of the quantum vortices. To see a quantum vortex, one 
must look for an isolated node of the wavefunction, i. e. \I/i/(2i) = sin(k x x) sm(k y y) + 
\l/ 2 (— x,?/)/(2i) = 0. Thus \& 2 ( — x, y)/(2i) must be real. Writing ^ 2 (~ x , y)/(2i) as 
R 2 exp(iS/h), we see that the nodes must lie on curves of constant phase S = integer xh/2. 
Now v = VS/m* would be the velocity if only ^ 2 (— x, y) were present, so constant S curves 
are orthogonal to v. Taking a region where the flow is to the right and downward, the 
orbits are straight lines at an angle —9 from the horizontal. Thus the constant S curves are 
also straight lines but at an angle 7r/2 — 9 upward. This is precisely the direction of the 
upward diagonal line of the standing wave cell. Picking a grid point of the standing wave, 
x = — (n/k x )n x , y = (ir/k y )n y , in the region we are considering now, 

V 2 (-x,y)/(2i) = - (1/ A) exp[i(k x + k y tan 9) ^n x ] x 

7T 71 

ep[i(—n y -—n x ta,n9)] (27) 

Now remembering that ep(iy) behaves as exp(ik y y) } 

* 2 (-x, y)/(2i) « -R 2 (-l)^y (28) 

which is real on these grid points. Therefore, the constant phase lines allowing quantum 
vortices are simply the upward diagonal lines of the standing wave cells. The line passing 
through the grid point (n x ,n y ) is 

7T z 71 " /in/ 71 " \ k x 7T 

y - TT n v = tan ( 9 _ ^(jT 7 ^ + x ) = TT^T Ux + X ) ' ( ' 

thy Z &> X l\jy K X 

Thus along the upward diagonal line, 

sm(k x x) sm(k y y) = (-l) n ^+ n « sm 2 (k x x) . (30) 

Therefore the condition to find a node is 

R 2 = sm 2 (k x x) . (31) 
9 



Of course, the solution will depend on the form of i?2- Once again, since ep(iy) ~ constant x 
exp(ik y y) in its no n- vanishing region, R 2 is a slowly varying function in the said region. Its 
average value along the diagonal line can be estimated from the transmission coefficient 
which yields l/|4cos#|, a value normally less than 1. Along the diagonal and within one 
cell, sm 2 (k x x) starts from zero at one end, monotonically increases to 1, then monotonically 
decreases to zero at the other end. Thus one normally sees two vortices per cell along the 
diagonal direction perpendicular to the average direction of motion of the particle. Fig. 3 
shows an array of such cells with two quantum vortices spinning in the opposite directions. 
The deviations of the vortex centers from the diagonal lines are small. Combined with Fig. 
1, one sees that the trajectories on average are mirror images of those in the transmitted 
region which resemble classical trajectories, but in detail they are much more complex. We 
have therefore found a complete justification for the ballistic trajectory, or the ray optical 
picture of propagation in the large wave number limit. Corrections to this picture will 
account for diffraction on both edges of the striped propagation pattern, and slow filling in 
of the shadow regions on both sides of the classical trajectory. This diffractive effect is seen 
in Fig. 1, and it is much bigger when the energy is lowered, as seen in Fig. 4. Notice that the 
average trajectories in the reflected region close to the classical trajectory rely heavily on the 
fact that the quantum vortices are situated according to the standing wave cell structure. 
The effect of a small misalignment of quantum vortices may be magnified significantly in the 
trajectories. This is because there is an equal number of hyperbolic unstable fixed points 
for each pair of vortices, and these can cause two orbits which are infinitesimally close in 
one place to follow quite different paths later on. This effect is also seen in Fig. 4. 
We would like to add the following additional points of discussion: 

• A natural question is, why do we not see turbulence despite the very complex flow 
pattern? The answer is that we are solving a steady flow problem. By definition, the 
velocity field does not vary with time. Since the stream lines do not cross each other, 
a single line can occupy very limited space. However, it would be interesting to see 
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what would happen if the steady state condition were lifted. 

• We have seen that as the energy is lowered, the trajectories start to deviate from the 
classical ones. An extreme is reached when the energy is so low that there is only one 
propagating mode. In this case, the wavefunction is well described by the single mode 
approximation || such that the transverse dependence is simply sm(iry/d) except in 
the region very close to the narrow constriction. Without the phase changing in the 
y-direction, the stream lines are parallel to the wire. As the stream lines approach 
the constriction, they will converge so that they may pass through the opening. After 
escaping, they then relax back and resume their motion parallel to the wire. This 
situation is common to all devices with straight leads, at low energy. 

• Since the stream lines do not cross each other, the bouncing ball trajectories do not 
really hit the walls. Rather, they turn to a more or less horizontal direction in the 
H-region( see Fig. 2), to avoid crossing. This is one difference from the strict classical 
prediction. 

• A larger opening width a will enlarge the H-region, and thus lead to stream lines more 
parallel to the wire as expected. 

• Previously we said that on the up-stream side, where not covered by \I/2 ( — x, y), there 
is no flow motion. This is an approximation. In reality, localized flow circulating in 
these regions can form quantum vortices, and possibly vortex clusters. This is a typical 
situation where stream lines close on themselves. 

• The problem solved in this paper can be taken as an odd symmetry solution for a wire 
of width 2d with a obstacle placed at its center, which is a variation of the ordinary 
double slits problem. Or as an odd-symmetry solution for a wire with a constriction 
standing out from each side of the channel. 

In summary, we have given a quantitative description of trajectories in a particular 
quantum wire transmission problem. The result closely resembles the classical bouncing ball 
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orbits in the large wave number limit. The major difference lies in the effect of quantum 
vortices which is responsible for the complex pattern in the up-stream. This simple model 
can serve as an explicit analytical example of the picture of ballistic transport in quantum 
waveguides, which in one or another form, clearly or vaguely, exists in our minds. 
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FIGURES 

FIG. 1. Stream lines participating in global transmission for a quantum wire with a simple 
constriction, a = OAd, E = 50.5 2 (7r/(f) 2 , ni = 40. The angle 6 reduced is 52.4°. 

FIG. 2. Formation of the bouncing ball pattern in the transmitted region due to the shift of 
the function ep(±iy) as x increases. 

FIG. 3. Current density plot in a portion of the up-stream propagating region. The grid 
standing wave pattern and associated quantum vortices are clearly seen. Horizontal and vertical 
lines are standing wave node lines. Inclined diagonal lines are the theoretical constant phase lines 
where VP(— x, y)/(2i) is real. Parameters are the same as Fig. 1. 

FIG. 4. Same as Fig. 1 but at a lower energy: E = 10.5 2 (7r/d) 2 , rij = 8. 
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